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Abstract 

We study early universe cosmologies derived from a scalar-tensor action 
containing cosmological constant terms and massless fields. The governing 
equations can be written as a dynamical system which contains no past or 
future asymptotic equilibrium states (i.e. no sources nor sinks). This leads to 
dynamics with very interesting mathematical behaviour such as the existence of 
heteroclinic cycles. The corresponding cosmologies have novel characteristics, 
including cyclical and bouncing behaviour possibly indicating chaos. We discuss 
the connection between these early universe cosmologies and those derived from 
the low-energy string effective action. 



1 Introduction 



In this paper we consider the qualitative dynamics of a class of spatially flat, scalar- 
tensor cosmological models derived from the action 
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where R is the Ricci curvature scalar of the space-time with metric g^, g = detg^ u , 
{A, Am} are constants and represent scalar fields. The dynamics of these 

cosmological models has very interesting mathematical properties. In particular, there 
are no asymptotically attracting equilibrium states in the phase space and this may 
lead to important physical consequences. 

The form of action (|l[) can be partially motivated from string theory, which is the 
most promising candidate for a unified theory of the fundamental interactions [|], 0. 
When Am vanishes, Eq. ([!]) represents the truncated effective four-dimensional action 
of the Neveu-Schwarz/Neveu-Schwarz (NS-NS) sector of the theory ||TJ. The scalar 
field, $, represents the dilaton field and the axion field, a, is the Poincare dual of the 
antisymmetric two-form potential. The constant, A, may be interpreted in terms of 
the central charge deficit of the string theory and can be negative. The evolution of the 
very early universe immediately below the string scale may have been determined by 
an effective action of this form. The dynamics of the spatially flat and homogeneous 
cosmologies in the case Am = was presented in Ref . . One of the main purposes of 
the present work is to determine the effects of introducing a non-trivial cosmological 
constant, Am, that does not couple directly to the dilaton field. Such a term represents 
a vacuum energy contribution to the energy-momentum tensor. A discussion of the 
spatially flat and homogeneous cosmologies in the case A = was presented in Ref. 

and a partial analysis of the dynamics with a non- vanishing axion field when both 
A and Am are non-zero was investigated in H. 



2 Analysis 



We assume that the metric corresponds to the spatially flat, Friedmann-Robertson- 
Walker (FRW) universe: ds 2 = —dt 2 + e 2a ^ dxidx 1 . Substituting this ansatz into the 
action ([[J) and integrating over the spatial variables then yields the reduced action 
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where a dot denotes differentiation with respect to cosmic time, t. The Friedmann 
constraint derived from Eq. (^) is given by 

' " * ' J - ,: ' A M e^ +3a = 0, (3) 
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in terms of the 'shifted' dilaton field, <p = $ — 3a. 

A generalization to the spatially flat, Bianchi type I cosmology may also be con- 
sidered. This effectively results in the introduction of two massless scalar fields into 
the reduced action (Q). These fields parametrize the shear of the models. Simi- 
lar degrees of freedom also arise when considering the toroidal compactification of 
higher-dimensional theories. Although we do not consider these extra fields in this 
paper, their overall contribution to the dynamics can be modelled by introducing a 
single modulus field, /3 2 = J2i$f, into the reduced action (0) 0], and their inclusion 
could be important in the discussion of chaotic behaviour. 

Zero central charge deficit 

We first consider the case A = |]. We assume that A M > 0, and employ the 
generalized Friedmann constraint equation (0) to eliminate the axion field from the 
system. The resulting field equations may then be simplified by introducing the new 
variables and time coordinate^: 

= _ -^m _d_ = ld_ = 1 {v+3a)/2 d_ m 

where a prime denotes differentiation with respect to 9 and ip = (p'. The Friedmann 
constraint yields 

l-x 2 -z>0, (5) 
from which it follows that the phase space is bounded with 

0<{x 2 ,z}<1. (6) 

The invariant set 1 — x 2 — z = corresponds to a trivial axion field. 

The cosmological field equations for the isotropic FRW model can now be ex- 
pressed in terms of the plane system: 



_ = ( X+y /3)[l-^- Z ] + l Z [ X -y^) ( 7 ) 



dx 
dO 

^ = 2z( [ [l-x 2 -z]-^(l-z-V3x)Y (8) 

The equilibrium points of this system 1,0), (1,0) and S + 

(-1/3V3, 16/27). The first two are saddles and S + is a repelling focus. The functional 
form of these solutions was presented and discussed in Ref. [|J and the phase portrait 
is given in Fig. 1. We note that the exact solutions corresponding to all of the 
equilibrium points are self-similar cosmological models 0. 

1 We assume that ^ > 0; the case ip < is related to a time-reversal of the system and the 
qualitative mathematical behaviour is similar (although the physical interpretation is quite different). 
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Figure 1: Phase portrait of the system corresponding to the isotropic FRW 

model with A M > and A = 0. We shall adopt the convention that large black dots 
represent sources (i.e., repellers), large grey-filled dots represent sinks (i.e., attrac- 
tors), and small black dots represent saddles. Note that in this phase space orbits are 
future asymptotic to a heteroclinic cycle. 



We see from Fig. 1 that the orbits are future asymptotic to a heteroclinic cycle. 
This is comprised of the two saddle equilibrium points Lf_\ and , and the single 
(boundary) orbits in the invariant sets z = (Am = 0) and 1 — x 2 — z = (a = 0) 
joining and W + y Hence, the orbits exhibit cyclical behaviour Q. For each 
'cycle', an orbit is quasi-stationary in the neighbourhood of the saddle point Lfy 
It then shadows the orbit in the invariant set z = as it moves rapidly towards the 
equilibrium point L* + y It settles into another quasi-stationary phase close to Lf + \ 
and eventually moves quickly back to shadowing the orbit in the invariant set 
1 — x 2 — z = 0. It is important to emphasize that the orbits move progressively closer 
towards the two saddles, Lf ± y after the completion of each cycle. Thus, the motion 
is not periodic and a given orbit spends more and more time in the neighbourhood 
of these equilibrium points. 

The physics behind the cyclical nature of these orbits is as follows. The sign of the 
variable x determines whether the universe is expanding or contracting. The value of 
this variable passes through zero during each cycle. This behaviour arises because the 
cosmological constant effectively resists the expansion of the universe, but the axion 
field has the opposite effect. Since the energy density of the latter scales as a 2 oc e~ 6 °, 
it is negligible when the spatial volume of the universe is large. Consequently, the 
cosmological constant forces the expanding universe to recollapse. However, the axion 
field inevitably becomes dominant and reverses this collapse, causing the universe to 
enter into a new expanding phase. The process is then repeated and the interplay 
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between the two opposing trends results in a universe that undergoes a series of 
bounces. 



Non-zero central charge deficit 

We now consider the case A 7^ 0. We again employ Eq. (^) to eliminate the a 2 
term from the field equations, and make the following definitions: 

_ V3a _ -2A _ A M e v+3a _ up d _ d 

We assume that A < (A M > 0) and we define £ 2 = ub 2 — 2 A. The generalized 
Friedmann constraint equation (|^) now yields < x 2 + z < 1 , so that all variables are 
bounded: < {x 2 , y, z, u 2 } < 1. From the definition of £, y is given by u 2 + y = 1. 
The resulting three-dimensional system therefore becomes: 

% = ^(l-i a -|jf)+zu(l-x 2 -± 2 ), (10) 

*= ( i _„>>(,» + !,)>„, (ii) 

= z [u (l - 2x 2 - z) + V3x\ . (12) 



The invariant sets x 2 + z = 1, z = 0, u 2 = 1 define the boundary of the phase 
space and it is important to note that the variable u is monotonically increasing. 
This ensures that there are no closed or recurrent orbits in the phase space. The 
equilibrium points of the system are all saddles: S ± (x,u,z = =pl/\/27, ±1, 16/27), 
L^n (±1, 1, 0) and LZn (±1, —1, 0). The points represent power-law cosmologies 
with ip > 0, where only the dilaton field is non-trivial, i.e., the axion field and 
cosmological constant terms are dynamically negligible. These solutions are termed 
'dilaton- vacuum' solutions and have an analytical form given by e a oc £±^3 and 
e* oc t~~ 1±v ^. The points are the corresponding solutions where ip < 0. The 
phase portrait is given in Fig. 2. 

In this case there are no sinks and no sources in the full three-dimensional phase 
space. Since the variable u (and hence (p) is monotonically increasing, solutions 
generically asymptote in both the past and future towards the invariant sets u — ±1. 
These both include an heteroclinic cycle and this implies that generically the solutions 
exhibit similar asymptotic behaviour at both early and late times to that discussed 
above. (For example, to the future the orbits in the three-dimensional phase space 
shadow the orbits in the two-dimensional invariant set A = 0.) The orbits interpolate 
between the dilaton-vacuum solutions corresponding to the saddle points in the 
past and the dilaton-vacuum solutions corresponding to the saddle points in the 
future. The effect on the dynamics of the cosmological constant, Am, is significant 
at both early and late times. The points S correspond to the equilibrium point 
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Figure 2: Phase diagram of the system lffi)-ftn\) for A < and Am > 0, where <p > 
is assumed. See caption to Figure 1. Grey lines represent typical trajectories found 
within the two-dimensional invariant sets, dashed black lines are those trajectories 
along the intersection of the invariant sets, and solid black lines are typical trajectories 
within the full three-dimensional phase space. 
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S + in Fig. 1, but unlike in the A = case in which S + is a repelling focus, in the 
three-dimensional phase space they are saddles and hence they do not play a primary 
role in the asymptotic behaviour. 

From a physical point of view, the cyclical behaviour arises due to the complex 
interplay between the axion field and the cosmological constant terms. The universe 
continues to undergo a succession of bounces between expanding and contracting 
phases due to the axion field and vacuum energy A M - However, the inclusion of the 
central charge deficit, A, causes (p to ultimately change sign. Thus, the asymptotic 
behaviour in the future is related to a time-reversal of the asymptotic behaviour in the 
past. We remark that the dilaton- vacuum solution, e a oc t^ 1 ^, corresponding to the 
point is inflationary over the range t < 0, because the expansion is accelerating. 
In this case, the accelaration is driven by the kinetic energy of the dilaton field. (For 
a recent review of the cosmological significance of these solutions see, e.g., Ref. [0]). 

Finally, we make some brief remarks on the case A > 0. We can define £ = ip 
and consider the subset ip > 0. Introducing normalized variables as before yields a 
three-dimensional, compact system of autonomous, ordinary differential equations. 
We have completed a full dynamical analysis of this system, but we only describe 
the main features here. There is a non- hyperbolic equilibrium point, C + , which can 
be shown to be a (global) source, since y is a monotonically decreasing function. 
This point represents a static universe, where the dilaton field is evolving linearly 
with time and the axion field and Am are dynamically negligible. There are also two 
saddle points, Lf^, which represent dilaton-vacuum solutions; these are analogues 
of the saddles that appear above. Again, there is also a saddle S + . We stress that 
there are no sinks in the phase space. Therefore, trajectories generically asymptote 
into the past towards C + , and then spiral away towards the heteroclinic cycle in the 
invariant set y — containing the saddle points Lt_\ and Lf + y The phase space is 
depicted in Fig. 3. 

3 Discussion 

The most important mathematical feature of the models we have considered is the 
cyclical behaviour that arises due to the existence of a heteroclinic cycle. This is 
of great physical significance, because it might be an indicator of possible chaotic 
behaviour. The solutions interpolate between different 'Kasner-like' dilaton-vacuum, 
power-law models, undertaking cycles between the saddles in the three-dimensional 
phase space. This is similar to the dynamical behaviour that occurs in spatially ho- 
mogeneous Bianchi cosmological models || |J. The question of chaos in anisotropic 
Bianchi type IX string cosmologies has been considered ||. It was shown that since 
the axion and dilaton fields behave collectively as a stiff perfect fluid, the system oscil- 
lates only a finite number of times. Consequently, there is no Mixmaster-type chaos in 
these models || |Hj] • This is to be expected since it is known that the admission of stiff 
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fluid matter causes chaos to cease |[L1]| . In contrast to the anisotropic Bianchi type IX 
cosmologies, however, the models described above contain cosmological terms, e.g., 
an effective dilaton potential and a dynamically important axion field. It is these 
intrinsically stringy effects that give rise to chaotic behaviour and this chaos has a 
different origin to the chaotic behaviour that arises in general relativistic models. On 
the other hand, there may be some connection with models that contain Yang-Mills 
fields. It is known that chaotic oscillations occur for such fields [fL3] and, moreover, it 



was shown in |TT] that the oscillations that are suppressed by a single massless scalar 
field can be restored by coupling an electromagnetic field to a Brans-Dicke type field. 
This model is related to a scalar field model with an exponential potential |13| and, 
consequently, is also related to string theory cosmological models [Q. 

There are a number of outstanding issues that need to be addressed regarding the 
possible existence of chaos in string cosmology. First, the chaotic behaviour depends 
crucially on the dimensionality of spacetime and on the product manifold structure of 
the extra dimensions 0. In particular, superstring theories are formulated in D = 10 
spacetime dimensions, while M-theory, with its low-energy supergravity limit, is an 
eleven-dimensional theory Second, the low energy effective action is only valid 
in the perturbative regime of weak coupling and small curvature. In general, it may 
be necessary to study chaos within the context of a full non-perturbative formulation 
of the theory, but at present such a formulation is unknown. Nevertheless, if chaotic 
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behaviour occurs at the level of the effective action, it is to be expected that similar 
behaviour should arise in the non-perturbative regime. Finally, there is the question 
of what will happen if inhomogeneities are introduced. Again, such effects will be 
most unlikely to lead to any suppression of chaotic behaviour and will perhaps make 
chaos even more predominant H. 

There are other questions which are important in early universe cosmology in gen- 
eral, and in string cosmology in particular. The questions of whether cosmological 
models can isotropize and/or inflate (and if they can inflate whether there is a graceful 
exit from inflation) are of great importance [|16||. The techniques utilized in this paper 



can easily be adapted to study the possible isotropization in more general spatially 
homogeneous but anisoptropic string cosmological models. Inflationary properties of 
simple string cosmologies have been discussed above. However, a chaotic cosmological 
regime might either be an alternative to inflation or, perhaps more importantly, could 
work in tandem with an inflationary mechanism [I7| to produce new interesting phys- 
ical phenomena. For example, a chaotic regime due to dissipative effects or chaotic 
mixing [|Ti|] could possibly be an alternative to inflation as a cause of homogenization 
and isotropization. This might alleviate the problems of initial conditions in inflation. 
This last point has been addressed in fl9 |, where it was suggested that there would be 
sufficient time for a compact, negatively-curved universe to homogenize since chaotic 
mixing smooths out primordial fluctuations in a pre-inflationary period. 
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